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We examine the two-point correlation functions of the fields exp(ia$) 
in the sine-Gordon theory at all values of the coupling constant (3. Using 
conformal perturbation theory, we write down explicit integral expressions 
for every order of the short distance expansion. Using a novel technique 
analagous to dimensional regularisation, we evaluate these integrals for the 
first few orders finding expressions in terms of generalised hypergeometric 
functions. From these derived expressions, we examine the limiting forms at 
the points where the sine-Gordon theory maps onto a doubled Ising and the 
Gross-Neveu SU(2) models. In this way we recover the known expansions of 
the spin and disorder fields about criticality in the Ising model and the well 
known Kosterlitz-Thouless flows in the Gross-Neveu SU(2) model. 



1 Introduction 



The most important generally unsolved problem in the subject of massive 
integrable quantum field theory in two dimensions is the computation of cor- 
relation functions. A completely general approach to the study of the large 
distance expansion is well known: one can insert a multiparticle resolution 
of the identity between fields and thus obtain an infinite integral represen- 
tation for the correlation function involving the form factors. Multiparticle 
form factors have been computed in a variety of models by Smirnov In 
principle, this form factor sum completely characterizes the correlation func- 
tion. However the complexity of the form factors makes this representation 
difficult to utilize. 

On the other hand, many quantum field theories have well defined short 
distance expansions. This short distance expansion is essentially indepen- 
dent of the large distance expansion discussed above; it is indeed extremely 
difficult to obtain one from the other. In this paper we develop the short dis- 
tance expansion of sine-Gordon correlation functions. This is also an infinite 
integral representation, though the integrals can be much more compactly 
written than the large distance expansion involving the form factors. It is 
our aim to demonstrate that this short distance expansion is tractable. 

The form of the two point correlation functions of the fields exp(ia$) in 
the sine-Gordon model at the free-fermion point (SGff) is well understood. 
Because these fields in the SGff can be mapped onto a doubling of the 
spin/disorder fields in the Ising model (Ityzkson and Zuber [|J), the results 
found by Wu et al.|| allow the expression of the correlators as solutions of 
a Painleve III non-linear differential equation. More recently, this result was 
extended in There, using a generalization of the techniques developed 
in ||, it was shown how correlators of the fields exp(ia$), < a < 1, are 
characterized by solutions to a non-linear equation, the sinh-Gordon equa- 
tion. The proof of this result lay in expressing the large distance expansion 
of the correlators as Fredholm determinants through the use of form factors, 
and then deriving differential equations for the Fredholm determinants. 

In this paper we present new results characterizing correlators of the 
above fields in conformal perturbation theory for values of the sine-Gordon 
coupling, (3, in the range 1 < (3 < \/2. (Here, (3 is related to the conventional 
sine-Gordon coupling constant as (3 = (3/ a/47t, the free fermion point occuring 
at f3 = 1.) Because of the U(l) charge symmetry of the conformal limit of 
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the sine-Gordon theory, we focus upon the following two correlators: 



Expanding in powers of the mass parameter, A, we are able to write integral 
expressions for every order of the short distance expansion. However, not 



correlators provided (3 < y/2. 

In principle it is possible to evaluate the integrals at every order. The sin- 
gularities appear as poles in the resulting functions: at lower orders, gamma 
functions, at higher orders, generalised hypergeometric functions. However, 
the results at higher orders quickly become intractable and not particularly il- 
luminating. Thus we restrict ourselves to evaluating the correlator, G(a, —a), 



to 0(A 2 ), and the correlator, G0/2,$/2), to 0(A 3 ). Indeed, G0/2,$/2) at 
0(A 3 J is already only finite for (3 2 > 3/2, one-half of the full range of interest, 
1< (3 2 < 2. 



The integral expressions in the perturbative expansion do not arise from 
a diagrammatic analysis. Thus it might seem necessary to explicitly subtract 
the contribution of the bubble diagrams. (Without such a subtraction, the 
terms in the expansion are formally divergent.) However, the technique, as 
first introduced by Dotsenko ||, we use to evaluate the terms renders them 
finite without explicit inclusion of the bubbles gous to di- 

mensional regularisation. Rather than analytically continuing the dimension 
of spacetime, we continue the parameter, (3. In this way not only are the 
integrals rendered finite, the bubble diagrams evaluate to zero. 

At particular points in the coupling, the sine-Gordon theory is known 
to map onto familar theories. At (3 = 1, sine-Gordon maps onto a doubled 
Ising model, and at (3 = y2, it maps onto the SU(2) Gross-Neveu model. 
We examine our perturbative expansion at these points. At both of these 
points the perturbative series become infinite and renormalisation is required. 
At (3 = 1 we introduce a novel renormalisation scheme to remove the IR 
divergences. At (3 = v2 we employ a more conventional scheme to take care 
of the UV divergences. Having then performed the necessary renormalisation, 
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we obtain the known expansions of correlators for the order /disorder fields 
in the Ising model as well as the known Kosterlitz-Thouless flows in the 
Gross-Neveu SU(2) model. This provides some degree of confidence that our 
derived expressions away from the free fermion point are correct. 



2 Overview of the Perturbative Expansion 

The action for the sine-Gordon theory in Euclidean space-time (our conven- 
tions are z — (t + ix)/2,z = (t — ix)/2 and d 2 z = —dtdx/2) is 

S = Soft + S Pert . = ~ j d 2 z(d z <bd- z § + 4A : cos0$) :) (2.1) 

where, again, the parameter (3 is related to the conventionally normalized 
coupling by (3 = -J=. As shown by Coleman JSJ, the UV divergences in 

the theory for values of (5 < \/2 (i.e. for values of (3 where the coupling A 
is not irrelavent), may be removed by the straight-forward normal ordering 
of the cos(/3<I>). In the process, the coupling constant A is multiplicatively 
renormalised. As no wavefunction renormalisation is necessary, the scaling 
dimensions of the fields remain unchanged by the introduction of the cosine 
perturbation - the dimensions are the same as those in the deep ultraviolet. 
Because the structure of the field theory is the same as its conformal limit, 
it is meaningful to perturb about the free theory - indeed, we may use the 
same labels for the fields in both the free and massive theories. For a general 
discussion of conformal perturbation theory see 

The coupling constant, A, in [O] can be directly related to the mass, m, of 
the asymptotic states, the solitons. Zamolodchikov |TT |, using results derived 



with thermodynamic Bethe ansatz, demonstrated that 

A = -A(/3)m 2 -^ (2.2) 



where [] 

A(/3) 



r(/3 2 /2) 



V5FT(l/(2 - /3 2 ; 



n 2-/3 2 



(2.3) 



r(i-/?72) L r(/? 2 /(4 - 2/? 2 ) J 

1 Some caution must be used in applying Zamolodchikov's result directly as he uses 
different conventions for both the action and the propagator. 
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We note that A = — m at (5 = 1. This relation ( |2.2| ) will prove to be useful 
in our examination of the Ising model. 
Consider the correlation functions 

G(a,a') = (0| : exp(io;$(z, z)) :: exp(ia'$(0)) : |0). (2.4) 

These correlation functions are highly non-trivial. Even at the free-fermion 
point, (3 = 1, the bosonisation relations, : exp(±$) := ^-t^, where the \I/'s 
are the spinor components in the massive Thirring model, imply the fields, 
exp(ia$), are not simply representable as fermions. 

The perturbative expansion of these correlators in a path integral formu- 
lation is 

G ^ = ~ jm^m < 2 - 5 ) 

J D$e~ ScFT( -^e* I d2w cos (/N™)) e ia*(z) e «*'<i>(o) 

/ D$e- s cFTm e $ J ' d2w CO80^(W)) 
00 1 / A \ n r 

E - ( - ) / # w i < cos(Mm)) ■ ■ ■ cos0^(w n ))e ia ^e ia '^ o ~> > C ft 

n=0 n - J 

— (-) / d2w i < cos(f3<5>(wi)) ■ ■ ■ cos([3$>{w n )) > C ft 
£±^A.y J # w . J2 < e «^*(^)... e ^*(^)e te *W e fa '*( )>ojr r 



n=0 

oo 



n=0 ' v 7 ffe=±l 



oo -j / \ \ n 



where we have expanded e~ Spert in powers of A and <>cft indicates the 
correlator is to be evaluated in the conformal limit. The advantage of ex- 
pressing G(a, a') in this form is that the correlators of vertex operators in 
the conformal limit, i.e. <>cft, are easily evaluated. The denominators of 
these expressions, representing the bubble diagrams in the theory, will prove 
to be important in proving the convergence of the terms in the perturbative 
expansion. 

To evaluate the resulting conformal correlators we need to specify the free 



boson propagator. With the action as given in 2.1 the propagator is 



($(2, z)$(0)) = - log(czz), (2.6) 
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where c is some arbitrary constant. As our convention we set c = 1. With 
this propagator the OPE of two vertex operators : e ia$< - z ' z ) :: e*^*^ :, is then 
equal to 

. e ia3>(z,z) .. e i/9*(0) ;= (\ z \ 2 y P : e i(«+0)*(*,*) . _| (2.7) 

Using this OPE, the terms of the form < e «*i*( z i) . . . e M "*' z "' >cft are then 
given by 

< e <o*»(«0 . . . >CFT= -Q ^ _ z .|2)^^' (2 g) 

That <>cft vanishes unless &i = is as consequence of the U(l) symme- 
try of the free massless theory. 

Given this U(l) symmetry, only certain values of a and a' lead to non- 
vanishing correlators, G(a, a'), in the massive theory in perturbation theory. 
We examine two in particular, G(a, —a) and G(/3/2, (3/2). The first has non- 
vanishing terms at all even orders in the perturbation, the second at all odd 
orders. One reason for choosing to examine these correlators in particular is 
their relation to the Ising spin and disorder field correlators, <cr(z, z)a(0) > 
and < fi(z, z)fi(0) >, respectively. At $ — 1, these correlators are related to 
the G's via 

2 . Mz,z) . MO). 1 /„, 1 1, „,1 1. 



< a(z, z)a(0) > 2 =< sin(^^) sin(-^) >= - (G(~, --) - G(~, - 

(2.9) 

and 

<M^^)M0)> 2 =<cos(^^)cos(^)>= ~ (g(± ~\) + G{\, -)) , 

(2.10) 

where we have used G(a, a') = G(—a, —a'). 

The perturbative expansion for G(a, —a) and G((3/2, (3/2), using [278| and 
the U(l) symmetry, simplifies to 



Iw- ) ( J / d 2 wi- ■ ■d 2 w 2n uj{(3,Wi)T(a,(3,w i ) 

G(a, -a) = "=o n "V H \ n J (2.11) 

00 1 ( \ \ (2n\ f 

5d U n J^---d 2 ^((3,w i ) 



n=0' 
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and 



00 1 / A \ 2n+1 r (2n + 1\ 

^ ~n\ [ 2tt ) J d2wi ' ' ' d2w2n+1 [ n _|_ 1 J 



l 2'2 J -1 ( \\ 2n (2n 



(2.12) 

where the functions a;, r, x, an d ^ are given by 



(A«; i )= II (K-^-lT 2 II (l^-^l 2 )- 42 , (2.13) 



l<i^j<n l<?<n<j<2n 



2i |2\a/3 



r(a,/3,^) = |^" 2 r n (\ Wl - z f\ Wl \-^ n (|^-z|- 2 |m 

l<i<n n<j<2n 

(2.14) 

x0,wi)= n (i^-^'iT 2 n (k-^iv 2 , (2.15) 

i<i^j< n l<i<n<j<2n+l 

2n + l>i^i>n 



ii>0,wi) = \z 2 \± n — ^i 2 ]^! 2 )^ n (| Wi _ z | 2 | w .| 2 ) " 2 . 

l<j<n n<i<2n+l 

(2.16) 

3 Convergence of Terms in Perturbative Ex- 
pansion 

In this section we examine the singularities that appear in the perturba- 
tive expansions of the correlators G(a, —a) and G([3/2, (3/2). We find that 
G(a, -a) at 0(A 2ri ) is UV finite for (3 2 < 2 and IR finite for (3 2 > 2 - 1/n, 
and that G0/2,$/2) at 0(A 2n+1 ) is again UV finite for (3 2 < 2 but IR finite 
for f3 2 > (2n + l)/(n + 1). We give a proof of this beginning with G(a, —a). 

To show UV finiteness, we consider the contribution to G(a, —a) at 
0(A 2n ). It takes the form 

d 2 Wl . . . d 2 W2n ( e ia $( z ) e ~ ia ®(0) Q-i^M . . . e 0*(w2n-l)g-0*(l«2n)^ 

— disconnected pieces. (3.1) 
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The leading order UV singularity occurs as W\ — > w 2 , — > 104, . . ., wm-\ — ► 
W2n- These singularities are governed by the operator product expansion 

. e ^»(») .. e -/3*fe) . = | x _ y |-2*» . e i$W»)-*(y)) . + finite pieces (3 2) 

= |a; - y|- 2/52 (7-/5 2 |x-2/| 2 : $,*(i/)fy*(y) :+■••), 

where in the second line pieces such as (x — y)d y $(y) and its antiholomorphic 
counterpart are discarded. Such pieces, when the OPE is inserted into the 
correlator, integrate out. 

The leading order singularity in this term does not come from the contri- 
bution of the identity operator in the above OPE; rather this term leads to a 
disconnected piece which is subsequently subtracted off. Thus the leading or- 
der singularity that contributes is \x — y| 2-2/3 . So the leading UV singularity 
at 0(A 2n ) in G(a, —a) has dimension 

d uv = n(2 - 2(3 2 + 2). (3.3) 

The n arises as there n approaches Wi — > Wi+i, the extra 2 from the integra- 
tion / d 2 (wi — w i+ i). Hence for p 1 < 2, i.e. d uv > 0, the terms are UV finite, 
as is expected from Coleman ||. 

The reader may worry that stronger divergences are introduced consider- 
ing the approaches w\ — * w 3 , w 5 — > w 7 , and so forth in addition to W\ — > w 2 , 
W3 — > W4, etc. But such additional singularities are governed by the OPE 

: d x <S>(x)d^(x) :: d v Q(y)dy$(y) := (3.4) 
\x - y\~ 4 (l + \x- y\ 2 : d y ${y)dy<&{y) + • 



Again the leading term in the OPE is disconnected. So its leading singular 
contribution is \x — y\~ 2 - But this is exactly cancelled by the integration 
/ d 2 (x — y) and so the overall UV singularity becomes no greater. 

When all the WiS are allowed to approach one another, successively sub- 
stituting the above OPEs into the correlator in pO leads one eventually to 
the correlator 

( e <«*(*) e -«*(°) : d w $(w)d^(w) :). (3.5) 

But this correlator, as can be easily checked, introduces no additional UV 
singularities as w approaches z or 0. Thus allowing in addition the w^s 
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to approach z or does not change the UV properties of the perturbative 
expansion. 

To demonstrate that G(a, —a) is IR finite at 0(X 2n ) for p 2 > 2 — 1/n, 
we recast each term in the perturbative expansion by making the conformal 
transformation 

Wi^l/wi. (3.6) 
The perturbative terms then appear as 

\w 1 \ 2 ^- 4 d 2 W 1 ■ ■ ■ \w 2 n\ 2 ^~ A d 2 W 2 n X 

/ e ia$(z) e -i«*(0) e i/3$(u;i) e -i/3$(w 2 ) . . . e i/3*(w2n-i)g-i/9*(w2n)\ 

— disconnected pieces. (3.7) 

The IR behaviour of this integral is determined by sending Wi — > for all 
i. The singular behaviour of the correlator in the above integral can be 
extracted by first sending w\ — > w 2 , w 3 — > w 4 , . . ., u>2n-i - ► ^2n, then 
W\ — > u> 3 , u> 5 — > W7, . . ., w 2n -3 — ► w 2n _i, and so on. Using the OPEs given 
in |!0 and ET4| and discarding, as before, the disconnected pieces associated 



with the identity, the dimension of the leading IR singularity is 



d 



IR 



2n (2{3 2 - 2)] + [2 - 2n(3 2 ] . (3.8) 



The latter term in dm is the contribution coming from the correlator, the 
former term from / \ wi\ 2 ^ ~ 4 d 2 Wi. It is important to not include the singu- 
larities that arise from the OPEs involving e ia *'°). These singularities are 
properly UV. Then for the term to be IR finite, we need dm > 0. Thus /3 2 
must be greater than 2 — 1/n. 

To check the UV finiteness of G(/3/2, (3/2), a similiar argument is used to 
the one above. At 0(A 2n+1 ) the term in the perturbative expansion appears 
as 

J d 2 w 1 ■ ■ ■ d 2 w 2n+ i x 

— disconnected pieces. (3.9) 

The most singular pieces of this expression come when w 2 — > w 3 , w 4 — > w 5 , 
. . ., w 2n —* w 2n+ i. These approaches are governed by the OPE in pO. Thus 
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the overall dimension of the leading UV singularity is 



d uv = n(4-2(3 2 ), (3.10) 

and so again for (3 2 < 2 the terms are UV finite. Additional approaches, say 
u>2 — ► W4, wq — > w$, etc., are govered by the OPE [O] and so, as before, do 
not change duv- If ai l the WiS, % > 1, are allowed to approach one another, 
say at w, the succession of OPEs leaves one to evaluate the correlator 

( e i$/2*(*) e 0mo) e -if>*{v>i) : d w $(w)d^(w) :). (3.11) 

But as w — > Wi, 0, or 2, no new UV singularities appear. 

To obtain the IR behaviour of G(/3/2, (3/2) for n > 0, we begin by treating 
the term at 0(A 2n+1 ) like the term at 0(A 2n ) in G(a, —a), contracting the last 
2n vertex operators e ±l @®^ in the correlator. However there is an additional 
contribution of 2f3 2 — 4 to dm in this case that comes from letting w — > w\ in 
3.11 . Thus the term is IR finite if (3 2 > (2n + l)/(n + 1). For n = it is easy 



to see that this formula still holds: the term at 0(A) is IR finite if p 2 > 1. 

Though we have identified the ranges of (3 where the perturbative expan- 
sion is divergent, we point out that the expressions we derive in the next 
section will only be divergent for individual values of (3: in general the ex- 
pressions we will derive will be meromorphic functions of f3. However it is 
not clear how to interpret the analytically continued expressions in regions 
of 13 where the original integral expressions are divergent. It is possible that 
that this continuation does not capture all the relavent physics. For example 
at f3 2 < 1, it may be the continued correlators do not reflect the presence of 
quantum breathers. 



4 Evaluation of Terms of Perturbative Series 

We begin by focusing on the integrals in the numerator of the perturbative 
expansion (i.e. we ignore the contribution of the bubble diagrams). Although 
these integrals are formally divergent, the techniques (as first developed by 
Dotsenko |Q) we employ to evaluate them render them finite in a manner 
somewhat analogous to dimensional regularisation. This handling of the in- 
finities suggests the technique takes into account the bubble diagrams. This 
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is born out by the fact that when the contribution of the bubbles are cal- 
culated explicitly using the same techniques, they turn out to be identically 
zero. 

The integrals we thus consider are of the form 

hn = J d 2 wi ■ ■ ■ d 2 w 2n u0, • • • , w 2n )r(a, f3, w 1: ■ • • , w 2n ), (4.1) 

and 

hn+i = J d 2 w 1 ■ ■ ■ d 2 w 2n +ix(P, wi, ■ ■ ■ , W2n+i)ip0, w ± , ■ ■ ■ , w 2n+ i). (4.2) 

These integrals are, in their full generality, unwieldly. Thus, in the end, we 
only explicitly evaluate I±, I 2 , and I 3 . Indeed, J 3 is already only finite for 
p 2 > 3/2, a part of the full range of interest. However, we will carry out the 
calculation for general n as far as seems reasonable. We begin with I 2n . 



4.1 Evaluation of / 



2n 



To make sense of I 2n we must first abandon light cone coordinates, writing 
it in the usual x-t coordinates. Doing so, and at the same time scaling out 
the z-dependence, we obtain 

hn = (-2) 2n \z\ 4n - 2a2+2 P 2n [dt 1 dx 1 ---dt 2n dx 2n x (4.3) 

n - tj? + ^ - x j) 2 Y n ((*< - hf + fa - x 3 f)~ p x 

l<i^j<n l<i<n<j<2n 

n (((^-i) 2 +^) 2 (^ + ^" XQ/3 

<i<n V 

n (^-if+x 2 y 2 {t 2 + x 2 ) 2 ) 



n<i<2n 



Now the branch points in hnS integrand of the Xj's lie on the imaginary 
axis. Thus the deformation of the x-contours pictured in Figure 1 is per- 
missible. This deformation corresponds to the following change of variables 
Xi — > —ie~ t2e Xi where e is some small positive parameter. Given that the 
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new contour 



Im(x) 



V 



Re(x) 



old contour 



Figure 1: A sketch of the deformation of the Xi contours. 



terms in the integral now take the form (U — tj) 2 — e l4e (xi — Xj) 2 , the fol- 
lowing changes of variables 



wf = U±Xi 



factorises these expressions leaving I 2n in the form 



'2n 



:-2) 



2n 



2n 



\4n-2a 2 +2f3 2 n 



X 



(4.4) 



J dwf ■ ■ ■ dw 2 h n dw 1 ■ ■ ■ dw 2n 6(wf, w£ n , e)9(w 1 , • • • , w 2n , -e) 
where 0(wi, • • • , w 2n , e) is given by 



6(w u ■ ■ •, w 2n ,e) = 



II (Wi-l- ie^) aP K - ie& i )- ttP x (4.5) 



a/3 



Li=l 



n . 2n 



j=i+l 
2n 



, *e(A,-A,)f II (w i -«; j -ie(A i -A J -))" 

j=n+l 



X 



2n 



a 2 



Y[ (wi-l-ieAi) al (wi-ieAi) al J| (wi - wj - ie(Aj - A.,-)) 

j=n+l j=i+l 



and Aj = iu+ — w, r . 

Having made this set of changes of variable, we are faced with a set of 
integrals that are represented by contours deformed around the branch points 
of algebraic functions. The signs of the coefficients of e in the integrand tell us 
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s w 

1 '2n 



"2n 



— • •: •: • — w. 

w " w " . . . . w " 1 1 

l 2 "3 »2n 

Figure 2: The marked points on each of the contours represent the branch 
points. We have not drawn the contours for w~,---,w~ n as they are ir- 
relavent to the argument being made. 

how to make the deformations. Focusing on the -u;~- contours, the coefficients 
of the relavent e's are functions of the variables wf. For various values 
of the wf's, the ^"-contours will enclose no poles, and the corresponding 
contribution to I 2n is 0. To find these values we consider a set wf ordered 
as follows: 

K > K > ■■■ > K n - 

We will show that if wf > 1 or wf < 0, the contribution to I 2n is zero. 

First suppose wf > 1. Then if we perform the w~ integrations in ascend- 
ing order of i, the set of contours for wj appear as in Figure 2. Because the 
u>~-contour can be closed at oo, the contribution to I 2n is 0. Only the branch 
points of w~ arising from the terms of the form (w~ — a + ief(wf, w^)) 1 are 
shown. The branch points of w~ arising from terms of the form (wf — a + 
ief(wf,w^))' y , i.e. terms where w^ appears multiplied by e, are along the 
imaginary axis of w~ and are taken to ±ioo as e — > 0. Thus they do not 
affect the deformation of the contour. 

Now suppose wf < 0. The set of contours for wf then appear as in 
Figure 3. Because the w 2n contour can be closed in the lower half plane, the 
contribution to I 2n is again 0. 

In addition to the branch points that arise directly from the terms (w~ — 
a+ief(wf , w^ )), branch points will also appear because of prior integrations. 
For example, after performing the first 2n-l ^"-integrations, I 2n is left in the 
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* 

M- 





* 

•— — w . 

1 '2n 



'2n-l 



•2n 



' w. w. 

■2 "3 

Figure 3: Pictured are the set of contours for w~ when w 2n < 0. The 
dashed lines indicate the contours may follow either path, depending on the 
specific values of the u>~'s. The *'s mark branch points arising from prior 
integrations, as explained in the text. 

form 

hn = J dwf ■ ■ -dw^dwzJiwf,- ■ ■ ,w£ n ,e) x (4.6) 

/ \ ~a.{3 / \ a/3 

(w 2 n ~ 1 + ieA 2n) [w 2n + ieA 2n ) g(w 2n -ie), 

where g is a generalised hypergeometric function with branch points at 
and 1. These branch points arise from integrating terms of the form 
(w^ — w 2n + ie(A k — Aj 2n )^) , k ^ i 2n - As the coefficent of e in this term 
is positive when = w~ n (as w£ — wf 2n > 0), g depends on w 2n — it as 
indicated. Hence the Wj 2n -contour flows underneath the branch points of g. 
The importance of all of this is, of course, that g's branch points do not 
interfere with the deformation of the contour that takes the contribution to 
0. 

Thus is has been shown that only for values of wf constrained to lie 
between and 1 is a non-zero contribution to I 2n made. As the are n! ways 
of ordering the tu+'s, there are n! cases to consider, each having the general 
form 

l> w^> ■■■> wt > 0. 

The contours for u>j~'s in this general case appear as in Figure 4. The branch 
points marked by the *'s now appear below the contours as we are performing 



•2n 
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* * 




Figure 4: Pictured are the positions of the set of contours for w i relative to 
the branch points when the wf are between and 1. 

the Wi j integrals in descending order of j. The placement of the w~. branch 
points beyond 1 is justified by the deformation of the contours that we now 
make as shown in Figure 5. As we see, the contours are now such that the 
values of the u>~'s are all greater than 1, thus justifying the placement of 
the w~ branch points beyond 1. However, the given ordering of these branch 
points should not be taken as fixed, but should be understood to depend 
upon the specific values of the u>~'s. For example, if we are at a point on 
the contour Ci such that w~ 2 > u>~ then for the contours Ck, k > 2 the 
branch point for w~ should always fall to the right of that of w~. Thus, for 
example, the ordering of the branch points for the contour C2n as pictured 
in Figure 5 assumes an integration region for the other 2n-l contours where 
w l> ,<••< w~. 

>'2n — l *i 

Now taking e — > 0, we thus are able to express I 2n as the following sum 
of integrals 

hn = H) 2 ^|" 2 -^ 2 ™ ]T / ((7(1), • • • , a(2n)) , (4.7) 
where S 2n is the permutation group of 2n objects and 

I (il, • • • , i2n) = Jl (il, • • • , i2n) -h (il, • • • , ^2n) , (4.8) 

where J\ and J 2 are given by 
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Figure 5: Pictured are the deformations of the contours employed when w 
through -U7+ are between and 1. 



■A (*!,•••, *2n) = tdw^W^il-W^ X (4.9) 
JO 

y o dw i2 u! i2 712 (i - w i2 ) 7 '2 (^ - w i2 y^ dw i3 • • • x 

fW2n-l 

/ ^ 2 „^ 2 „ (1 - ^J^" K - ^ 2 „) 7<li2 " • • • 

v 

J 2 (ii,---,^n) = f dwi.w'^il-w^x (4.10) 

/ dw i2 ^ T<2 (l - »i 2 ) 7t2 Ki - W i2 ) 7 'i l 2 / G^ 3 • • • X 

[ dw i2n w~l l2n (i - w i2 y^ (w h - w i2n y^n ... x 

J Cm 

( on. _ on. YlHn-l^n 

and we have defined the 7's as follows 

r = l a(3/2 l - l ~ n ■ (411) 

h \ -a/3/2 n < i < 2n ' v ' ; 

/ P 2 1 < i, j < n or n < i,j < 2n 

7*3 ~ ) _n2 l<i<n l<j<n ■ V±- 11 ) 

I l J n<j<2n UI n<i<2n 
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-l Try: 
e 1 



w. 



e 1 



o V — W x — 12 

* * e n 2 w j e W1 2 'hh 

Figure 6: Shown are the phase assignments for the integral J^iHi^)- 



In summing over the Ji^'s we rearrange the signs of the terms involving 
powers of Wi — Wj. But because we always do this in pairs, one for J\ and one 
for J 2 , there are no phases between the terms in the sum J2aeS n - We note 
that . . . , ?2n) is invariant under various permutations of the i/s. Any 
permutation not mixing values of the i/s between 1 through n with those 
n+1 through 2n leaves . . . , i 2n ) unchanged. 

Now J\ is no more than a generalised Euler integral which we evaluate in 
Appendix A. The form is sufficiently complicated for the general case that 
there is little point in writing it down here. Indeed, ^]9| offers a much compact 
expression for J\. on the other hand, is as of yet still a formal expression. 
The contours Ck enclose multiple branch points and the corresponding cuts 
introduce phases. These phases depend upon which section of the integration 
region we are in, as explained earlier. These sections can be labelled by the 
specific order of the u^~'s in the section. As there are n! ways of ordering 
the w~'s, J 2 is a sum of n! terms. Now it is possible to develop a general 
notation for the assignment of the phases in each of the n! cases. However, 
this does not bring us any farther in actually evaluating I 2n . Thus it is at 
this point that we specialize to specific values of n, in particular n = 1. This 
case, together with the evaluation of ^3, brings out all the subtleties inherit 
in the assignment of phases. 
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4.2 Evaluation of I 2 

Assigning phases is particularly straightforward for the case of I<i- From 
Figure 5, n=l yields phases as pictured in Figure 6. ^(ii,^) then has the 
following form: 



where we have introduced the notation sia) = sin(7ra). In writing this ex- 
pression for ^(ii, %2) down we have assumed that the portions of the contours 
that circle about the branch point contribute nothing. This assumption can 
always be made good by continuing the exponents of the terms in the inte- 
grands to points where these portions do contribute nothing, evaluating the 
integral, and then continuing back to the desired value of the exponents in 
the final expression. It is this continuation of the exponents that makes the 
method akin to dimensional regularisation. But instead of continuing the 
dimension of space-time, we continue the the parameter (3, i.e. the scaling 
dimension of the vertex operators appearing in the action. 

^2(^1,^2) is now in a form where it may be directly evaluated. Making 
the change of variables, Wi — > l/u>j, the above may be written as 

•Mil, k) = -4 [s(7*i)s(7«)-K'(i2, h) + s(7*i)s(7ia + 7«i 2 )- f ^i, ^)] , (4.14) 
where K(i,j) is defined to be 



K(i,j) in Appendix B (even though the K's are no more than Euler integrals, 
Appendix A is not sufficient to evaluate the K's) with the result 




(4.13) 





K(i u i 2 ) = B(i + lh J 2 -i)r(i-/f)r(/3 2 - i)(- 7i2 )(/? 2 



1) 




;2 -2 



(4.16) 
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where 3F2 is a standard hypergeometric function. 

Using Appendix A to directly evaluate Ji(ii,i 2 ) we find 

Ji(ii,i 2 ) = J o dw^w'^il -w^h x 

[ n dw l2 w- % * (1 - w i2 p* { Wil - w i2 )™ (4.17) 

= £(l + 7il ,2-/3 2 )5(l-/3 2 ,l- 7i2 )x 

3 F 2 (- 7 , 2 , 2 - /3 2 , 1 - 7i2 , 3 - /3 2 + 7n, 2 - /3 2 - 7l2 , 1). 

So -Z2 has the form 



-A\z\^ 2a2 - 2 ^ 2 [^(l, 2) (s 2 (f )^(2, 1) + s(f )s(/3 2 + f )X(1, 2 

+j x (2, i) ( s 2 (f )K(i, 2) + S (f ) S (f - hm i)) 

(4.18) 



4.3 Evaluation of / 



2n+l 



We now move on the I 2n+ \. All of the same techniques that were used to 
analyze I 2n apply here. We are thus able to immediately write I 2n+ i as 
follows: 

h n+l = {-if n+ M 2{2n+l) -^ 2n+i) E /(<r(l),...,<7(2n + l)), (4.19) 

where I = (h, ■■■ , i 2n +i) = -h (h, hn+i) J2 (h, • • • , hn+i) , and J x and J 2 
are given by 

= [ 1 dw ll wi n (l-w H ) 5 ^ x (4.20) 

v 



y *3 



rfw i2 wf; 2 (1 - u>i 2 )^ 2 - u>* 2 ) 5l112 / "dwx ■ ■ x 



- 12 -jo 

f W2n S- A A 

/ dw^y)^ 1 (1 - ^ 2n+ i) l2 " +1 (w tl il<2n+1 
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-i Jt8i 
e 1 



1 i7l8, 
e 1 



Figure 7: Shown are the phase assignments for the integral ^(l)- 

h (il, • • • , «2n) = / dw^Wi* 1 (1 - lUiJ^i x (4.21) 

f dw i2 wf; 2 (1 - w i2 ) 5i * (w h - w i2 ) 5i ^ f dw h • • • x 

J Co. J C-x 



J C2n + l 



s5i 



r 

...(in- — 111- P'2n*2n+1 

and we have defined the <5's as follows 

^ = {7f/ 2 !<<<» + ! (4 . 22) 
\ /5 2 n + K i < 2n + 1 V ; 

and 

Oij — i _o2 l<i<n+l l<i<"+l 
I ' J n+Kj<2n+l n+Ki<2n+l 



(4.23) 

Again J 2 is a formal object because as of yet phases have not yet been 
assigned to the various portions of the contours C^. And again there is no 
advantage of performing this assignment in general and so we specialize to 
the cases of h and J 3 . 

4.4 Evaluation of I\ and /3 

Now Ii is particularly simple. The phase assignment is given in Figure 7. 
Thus J2(l) equals 

J 2 (l) = 2is(5 1 ) J™ dw^iw! - l) Sl (4.24) 

= 2is(£i) f dw lW i 2 - 25l (l -W!) 51 . 
Jo 
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Now Ji(l) and J2(l) are no more than the integral representations of the 
beta-function. Thus 

Ji(l) = + 1 + 50, (4-25) 
J 2 (l) = 2is(5i)B(-l-2(Ji,l + (Ji), (4.26) 

and so because 5i = — /3 2 /2 

/i = -2|^| 2 ^ 2 / 2 S (^/2)B(1 - y , 1 - y)B03 2 -1,1- f ). (4.27) 

We mention that the integral I\ is one that also arises in calculating the 
4-point Virasoro-Sharpiro formula. The methods used in this case for evalu- 
ating this integral give the same answer as above. 

The evaluation of J 3 is more complicated because the assignment of phases 
is more complicated. The assignment differs depending upon the ordering of 
tUjj with Wi 2 . The phase assignment for the two possible orderings is pictured 
in Figure 8. Defining K (11,12,13) by 



K(i ± ,i 2 ,i 3 ) = J dw^w^iw^-l) 8 ^ x (4.28) 
J 11 dw i2 w i2 s ^{w i2 - l) Si2 (w h - w i2 Y 1 ^ x 

dw i3 w i3 Si 3(w i3 - l) Si ^(w h - Wi 3 ) 5 ^3(w i2 - w i3 ) 5 ^, 

which after a change of variables Wi — > l/wi can be written as 

K(i u i 2 ,i3) = f 1 dw i3 w i3 - 2 - 25i s- 5 ^- 5 hh(i- Wi )* 3 (4.29) 
Jo 

[ 13 dw i2 w i2 - 2 - 2Si 2- Si ™- Si *n (i _ Wi2 ) s *2( w - Wi2 ) s ^2 
Jo 

P 2 dw h w h ' 2 ~ 25 ^ s ^- s ^(l - w h ) s ^ 
Jo 



W i3 -WiJ^n x 



(w i2 - w^) 6 *™. 

allows us to write J 2 (i\, ii, h) as 
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'1 »2 



-i 7lS: 

e l 

c 



o V j w. 

v_y ,j 

1 i7i& 

* * e 1 
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• — ~ • ,- — ■ — ; V, 

Figure 8: Shown are the phase assignments for the integral hih, «2, h)- 
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J 2 (ii,M 3 ) = (2i) 3 [ S (5 il ) S (^ 2 ) S (5 i3 )X(i 1 ,i 2 ,i 3 )+ (4.30) 

s(S il )s(8 i2 )s(8 i3 + 5 i[ii2 )K(i u i 3 , i 2 ) + 
s(S h )s(8 i2 )s(8 i3 + S i3i2 + 5 i3il )K(i 3 , i u i 2 ) + 
s(5 h )s(5 i2 + 5 ili2 )s(5i 3 )K(i 2 ,ii,i3) + 
s(5 h )s(5 i2 + 5 ili2 )s(5 i3 + 8 i3il )K(i 2 , i 3 , k) + 
s(8 h )s(5 i2 + 8 hi2 )s(8 i3 + 8 i3h + 8 i3 i 2 )K(i 3 , i 2 , h)} . 

Because K(i±,i 2 , i 3 ) is invariant under permutations of the i's involving 1 and 
2, it can easily be shown from the above expression that 

J 2 (l,2,3) = 0; (4.31) 
J 2 (l,3,2) = (2i)V(/3 2 /2)( S (3^/2)+ S 2 /2))X(l,2,3); (4.32) 
J 2 (3,l,2) = (2 ? ) 3 S (/3 2 /2) S (3/3 2 /2)( S (3/? 2 /2) + S (/? 2 /2))^(l,2,3) + 

s 2 2 /2)(s0 2 /2) + s(3/3 2 /2))X(l, 3, 2). (4.33) 

J 2 (l, 2, 3) vanishes indentically because of a cancellation of the phases. Using 
the fact that both J\ and J 2 possess the same symmetry as K, I 3 then takes 
the form 

h = 16M 6 - 7 ^/ 2 [K(l, 2, 3) (s 2 (/3 2 /2)( S (3/3 2 /2) + S (^ 2 /2)) ^(1, 3, 2) + 
S (/3 2 /2) s (3/3 2 /2)( S (3/3 2 /2) + s (/3 2 /2))J 1 (3, 1,2)) + 
K(l, 3, 2) ( S 2 (/3 2 /2)( S (3/3 2 /2) + s{fr/2)) J : (3, 1, 2) 



(4.34) 



Now all that remains is to evaluate the Ji's and the K's. 
Using Appendix A we find 



K(l, 2, 3) = B(l + fi 2 /2, -3 + 2/3 2 )B(l - /3 2 , -2 + 3/3 2 )B(l + /3 2 , /3 2 - 1) 

(/? 2 /2) fcl (/3 2 /2) fc2 (/3 2 ) fe3 (-3 + 2/? 2 ) 

fcl+fc2 



X 



oo 



E 



k f£ k3 h\k 2 \k 3 \ (_ 2 + 5/5 2 /2) fel+fe2 

x (~ 2 + 3 / 3 2 )fc 1 +fc 2 +fc 3 (-1 + /3 2 )fc 2 +fc 3 . (435) 
(-1 + 2(3 2 ) kl+k2+ k 3 (2(3 2 )k 2 +k [i 
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K(l, 3, 2) = B(l - (3 2 /2, -3 + 2(3 2 )T(1 - /3 2 )r(-2 + (3 2 )T(1 - (3 2 )T((3 2 - 1) 
" (-py 2 ) kl 2 /2) k2 (-p 2 )k 3 (S + 2p 2 ) kl+k2 



x 



E 



klk2k3 hlW. (-2 + 3p 2 /2) kl+k2 

(~ 2 + P )k 1 +k 2 +k 3 (-l+/3)fc 2 +fc 3 , s 

r(-l + A* + fc 2 + fc 3 ) T(A; 2 + fc 3 ) ' 1 ' J 
Ji(l, 3, 2) = B(l - /3 2 /2, 3 - 3/3 2 /2)B(l - /3 2 , 2 - /3 2 )B(1 - /3 2 , 1 - /3 2 /2) 

)ki+k 2 

klk2k3 h\k 2 \k 3 \ ~ (4 - 2p 2 ) kl+k2 

x (2-/3 ) fcl+ fc 2+ fc 3 (1 - /2)fc 2 +fc 3 _ (437) 
(3 - 2(3 2 ) k l+k 2 +k 3 (2 - 3/3 2 /2) 
Ji(3, 1, 2) = B(l + /3 2 /2, 3 - 3/3 2 /2)B(l - /3 2 , 2)B(1 + /3 2 , 1 - /3 2 /2) 
x ~ 2 /2) kl 2 /2) k2 2 )k 3 (3 - 3/3 2 /2) fcl+fca 
h\k 2 \k 3 \ (4-/3 2 ) fel+fe2 

x (2) fcl+fc2+fc3 (1 - [3 2 /2) k2+k[i (438) 
(3 - P 2 ) 

ki+k 2 +k 3 

(2 + [3 2 /2) k2+k3 

In evaluating these integrals the same sort of techniques that went into eval- 
uating K{i\,i 2 ) (i.e. handling the appearance of a pole multiplying a zero - 
see Appendix B) were employed. This completes the evaluation of I 3 . 

4.5 Evaluation of the Bubbles 

We now examine the contribution of the bubble diagrams (i.e. the integrals 
in the denominator of the perturbative expansion). We will show that with 
our techniques for evaluating the integrals, the bubbles vanish identically. 
The integrals that we must evaluate are 

/&= fd 2 w 1 ---d 2 w 2n n (K-^iT 2 n (k-^-iv 2 - (4.39) 

J i<i^j< n l<i<n<j<2n 

2n>i^j>n 



x g (-p 2 /2) kl (P 2 /2) k2 (-P 2 )k 3 (3-3/3 2 /2), 



Performing the change of variables i — ti ± Xi as before, 

omes 

l 2n = (~i) 2n Jdwtdw-r](wt,---,wt n ,e)r](w^,---,w 2n , -e), 



I 2n becomes 
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1 



w. 



'1 



■2n-l 



l 2n 



Figure 9: Pictured are the set of contours for w~ arising from the integration 
of the integral I^ n . In contrast to the situation for I 2n and hn+i, there are 
no branch points at and 1. 

where rj takes the from 

2n n 

rj{ Wl , • • • , w 2n , e) = J] II K - ^ - ie(Ai - A,)) ±/52 , (4.40) 

j=l j=i+l 

where the sign of the exponents deponds upon the particular values of % and 
j. Again taking the ordering of the w^s to be > ■ ■ ■ wf 2n1 the deformation 
of the w~ contours about the branch points the appear as in Figure 9. 

Because there are no branch points at and 1, the ^"-contour (or, in- 
deed, the w~ n -contour) can always be closed at oo. Thus is identically 
0. In Figure 9 a specific ordering of the branch points has been as- 
sumed. However, changing this ordering does not change the argument just 
made. The vanishing of the bubbles is supported by our analysis of the Ising 
spin correlators. Because our results match previous calculations, it is likely 
that we have accurately taken into account the perturbative structure of the 
vacuum. 
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4.6 Summary of Results 

At this point we gather previous results to provide expressions for G(a, —a) 
and G0/20/2): 



G(a, —a) 



\-2a 2 



(4.41) 



-2a 2 - 



X 



|4-2o 2 -2/3 2 



X 



JxCl, 2) (s 2 (f )K(2, 1) + s(f )s0 2 + f)K(l, 2)) + 



^1(2,1; 



S 2 (f)iT(l,2) + S 



A 

2^r 
A 

7T 



(4.42) 



2-3/3 2 /2 



+ 



A 



7T 



s(/3 2 /2)B(l - /3 2 /2, 1 - (3 2 /2)B0 2 - 1, 1 - /3 2 /2) 

z |6-7,3 2 /2 x 

2, 3) ( S 2 (/3 2 /2)( S (3^ 2 /2) + S (/5 2 /2)) ^(1, 3, 2) + 
S (/3 2 /2) S (3/3 2 /2)( s (3/3 2 /2) + S (/3 2 /2)) ^(3, 1, 2)) + 
K(l, 3, 2) ( S 2 (/3 2 /2)( S (3/3 2 /2) + S (^ 2 /2)) ^(3, 1, 2 



The expressions for J(i,j) and K(i,j) are given in |4.16| and [4. 17| respectively. 
The expressions for J(i, j, k) and K(i,j, k) are given in [4.35| through [4.38|. 



5 Sine-Gordon as Ising and SU(2) Gross-Neveu 

In this section we examine our expressions in the limits in which the sine- 
Gordon theory maps onto a doubled Ising model ((3=1) and the Gross-Neveu 
SU(2) model = y/2). We first examine the doubled Ising model. 



26 



5.1 Doubled Ising Model 

As demonstrated in section 3, as (5 approaches 1, the perturbative expansion 
develops IR singularities. In our expressions, these singularities manifest 
themselves as poles in gamma functions. These poles are indicative of the 
logs that appear in the Ising spin-correlators in the scaling limit (see, for 
example, ||). Specifically, a n-th order pole term indicates a n-th order log. 
The logs are obtained by expanding the powers of \z\ multiplying the pole 
terms. For example, writing (3 2 = 1 + e with e small, a typical situation will 
see the following expansion: 

uM 2 i 

log |«| H . (5.1) 



(3 2 - 1 e 

After this expansion is made, the divergent pieces (i.e. the 1/e term) are 
subtracted away. This can be understood in analogy with a cut-off method 
of regularisation. Throwing away the divergent terms amounts to no more 
than a specific choice of the cutoff. This procedure differs from that used 
by Dotsenko §| where the powers of \z\ are not expanded and logs are sub- 
stituted directly for the 1/e terms. This procedure is equivalent to ours in 
situations involving first order logs. However it fails to produce the correct 
coefficients of higher order logs that arise from overlapping divergences. 

To test this regularisation technique, we examine a correlator simpler 
than the Ising correlator, the free-fermion correlator. As such we will need 
to review some aspects of the connection between the free-fermion and the 
sine-Gordon theory at j3 — 1. The free fermion action in Euclidean space- 
time reads 



S — — / dxdt 

4:71 



ip^d z ip + + ip-dzip+ + im(ip-ip+ — , (5.2) 



where we have introduced the Dirac spinors ty± = (^) which satisfy \& + = 

and employed an appropriate choice of gamma matrices Q. The equations 
of motion for these fields are 

d z ip± = irmp±, 

(5.3) 

dz?p± = —imip±. 
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The correlators we will examine are 

$_(z,^ + (0)) and (^_(z,z)V+(0)). (5.4) 

Knowing the first of these correlators, the equations of motion allow us to 
fix the second. With this in mind, the on-shell mode expansions are: 

V>± = ±V^ r — e 6 ' 2 (c*(e)e^ 6 > x - d ± (9)e ±i ^> x ) 

(5.5) 

$ ± = -zV^ [°° — e- e / 2 (c^{e)e^ e > x + d h {6)e ±i P ( - e > x ). 

In writing these expansions we have introduce the rapidity variable 9 by 
which on-shell Euclidean momentum-energy is parameterized, 

E = imcosh(9) and p = msinh(6'), (5-6) 

so that — E 2 — p 2 = m 2 and ip{9) ■ x = m(ze e + ze~ e ). The creation (+) and 
destruction (-) operators obey the following commutation relations 

{d~(9),d + (9')} = {c-(6),c + (6')} = A<k 2 5{9 -9'). (5.7) 

Using these mode expansions together with the equations of motion, it is 
then straightforward to calculate the two correlators 

(ip-(z,z)ip + {0)} = -2imK (mr), (5.8) 

(V>-(z,z)V>+(0)) = 2my5K 1 (mr). (5.9) 

where the K's are the standard modified Bessel functions. 

The relavence of these correlators is found in the demonstration of Man- 
dlestam ]TIJ that the spinor components can be expressed in terms of the 
quasi-chiral components of the sine-Gordon field, <p L and <p R : 

ip± = exp(±i(f) L ), 

(5.10) 
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The fields cj) L and <p R are termed quasi-chiral because in the massless limit, 

d- z <p L = d z <f) R = 0. (5.11) 

In this limit (p L is the left mover and <fi R the right mover. 

With these conventions we can make the following identification 

(e^ (f V^ (0) ) = -2imK (mr). (5.12) 

At first order the l.h.s. of this equation reduces to 



F0=1) = — I d'wie^^e^^e-^'^cF, 



2tt ' " 

A f , 9 1 1 , 

= -i— / d 2 w- --. 5.13 

The i arises from taking apart e**, i.e. e®( w ' w ^ = — Iq-^ i w ) e - l <t> R ( w ) _ 

This integral has a logarithmic divergence. To regulate it we continue the 
scaling dimensions of the operators (via continuing p) as follows: 

F 0) = — [ tf w ( e ifo*(*) e ifi*(o) e -0*{«>,*)) 

= _,A f Sw — 7t- — (5.14) 

2tt J (w-zy 2 wP 2 V ; 

Notice that the U(l) charge of the vertex operators still sums to zero. Writing 
this integral in terms of Euclidean x-t coordinates, we easily evaluate it to 
be 

F0) = — kr-^ g( o 2 ~ o 2 f 2 2) B(2/3 2 - 2, 2 - n (5-15) 
7T 2 - 2(3 2 

We now intend to take /3 2 to 1 via /3 2 = 1 + e. A then has scaling 
dimensions of 1 — e. So we may express A as A = — m/i _e where fi is a mass 
scale fixed to be 

fi~ e = (l-e7 + 0(e 2 ))m~ £ (5.16) 

by Zamolodchikov's A — m relationship (see ^.2| ). Taking e — ► in F(p) then 
gives us 

F(l + e) = --m^ e \zfx\- 2e (5.17) 
e 

= -imtf Q - 2 log(m|z|) - 2 7 ^) (5.18) 
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As indicated before we drop the 1/e term. The above expression then reduces 
to 

(TfiT \ 
log(— )+7j (5.19) 

as \z\ = r/2. In taking the limit e -> we have not expanded out the // 
term so as to preserve the scaling dimension of F{p). Expanding out K (mr) 
to first order in m shows the two sides of |5.12| to be in agreement. We thus 



have gained some confidence that our methods of evaluating the perturbative 
terms are correct. 

To illustrate the connection between our method of regularisation and a 
cut-off method of regularisation, considering regulating F(l) with a cut-off 
\w\ < R: 

F(1) = -<A ffajJ—l. (5.20) 

Z7T J [W — Z)W 

Doing the integral we find 

R 

F(l) = 2i\ log(— ) + const. (5.21) 

\z\ 

Now R can be written as R = abm~ l where a and b are dimensionless 
constants and a is to be taken to oo (as the cutoff is removed). So as A = — m 

F(l) = 2zmlog(crV 1 m|z|) + const. (5.22) 

To obtain a match between F(l) and the known form of (ip-(z, z)ip+{0)) 
log(a) must be thrown away and a specific choice for b must be made. The 
throwing away of log(a) corresponds to our discarding 1/e. However in our 
method, to our advantage, there is no need to fix the arbitrary constant, b. 

The fermion correlators also provide an opportunity to test whether we 
are accurately taking into account the vacuum structure of the theory. The 
correlator (ip_(z, z)tp + (0)) has its first non-trivial contribution at second or- 
der in A, the first order at which the bubbles contribute. Bosonizing this 
corrrelator and expanding to second order we obtain 

2mW-Ki(mr) = ($-(z, z)i>+(0)) (5.23) 



z 



- + -Lk(1) + 0(\ 4 " 
z At: 2 v 
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where 



A 2 

K(l) = — / d 2 w x d 2 w 2 w x w^ 1 (w x - zY 1 (w 2 - z)\w x - w 2 \~ 2 . (5.24) 



We note that the zeroth order term of this correlator agrees with expansion 
of Ki. 

The integral K in |5.24| is divergent. (In fact it is quadratically divergent. 
If the bubbles were explicitly included we would only be facing a logarithmic 
divergence. That our evaluation of K still leads only to log's indicates we 
have correctly handled the bubbles.) Regulating K{1) as before by continuing 
j3 we obtain 

\2| 14— 2/3 2 

K0) = / d 2 w x d 2 w 2 w~fw 2 p \w x - l)-?\w 2 - lf\w x - w 2 \^\ 



K x 0) = [dwj dw 2 {w x - w 2 )-P = \ ^ } (5.27) 



(5.25) 

Using the techniques described in section 4, K(j3) is easily evaluated to be 
K0) = -AX 2 \z\ A - 2 ^z-^ x (5.26) 
K x 0) (s 2 2 )K20) + (s0 2 )sW 2 ) + s 2 2 )) K 2 (-P) 
where K\{p) and K 2 (±f3) equal 

r«i „ B(l -/3 2 ,1) 

h 'Jo " ' " 2-0 

K 2 (±P) = f'dw.wt^il-wt)^ 2 r dw 2 w'f- 2 {l-w 2 )^\w x -w 2 )-^ 
Jo Jo 

= p 2 (i - p 2 )r(P 2 -1,1- P 2 MP 2 - 1, i t /? 2 ) 

3 F 2 (l T /3 2 ,/3 2 ,/3 2 -l,2,/3 2 T /3 2 ,l) 

= /3 2 (i-/3 2 )r 2 (/3 2 -i)r(i-/? 2 )r(2-/3 2 ) 

3 F 2 (l±/3 2 ,/3 2 -l,2-/3 2 ,2,l,l), (5.28) 

where in the last line we have used the analytic continuation formula in 
appendix C. Now writing f3 2 = 1 + e and taking e — ► we find 

#i09) = -^(l + e + 0(e 2 )) (5.29) 
K 2 0) = l( e + 0(e 2 )) (5.30) 
K 2 (-P) = ^(l + e + 0(e 2 )). (5.31) 
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Some of the details of this calculation may be found in appendix D. Com- 
bining these results we find K(l + e) to be 

K(l + e) = -47rV^lV| 2 ~ 2/§ V^(l + e) 

= 4vr 2 M V" 2 (-1 + 2 (log(^) + t - I)) (5.32) 

where we have used /z~ e = (1 — e^j + 0(e 2 )) m _e and \z\ = r/2. We again pre- 
serve the correct scaling dimensions of K{(3) by not expanding z~^ ' . Instead 
we expand only |z| 4_2/32 , the z-dependence scaled out from the integral. It 
is this piece, coming from the integral, that we expect to generate the log's. 
Dropping the infinite piece in K(l + e), we find for (i/)-(z, z)ip + (0)) 

(V-M)V>+(0)> = - + -\mz\ 2 [\og(— ) +7 --j+0(m 4 ) 



/z ( 1 mr / ,mr . 1\ , . 



(5.33) 



Expanding out Ki(mr) } we find the two sides of |5.23j agree. Thus the bubble 
diagrams, even though not explicitly included, are taken into account. 

Having demonstrated that our regularisation techniques work on fermion 
correlators, we now apply the same methods on the Ising spin correlator with 
some confidence. We recall the spin correlator is given by 

^V(0)> 2 = \ -\) - C?(i -)) . (5.34) 

From Wu et al. , the expansion in the scaling limit of the spin correlator 
is 

(a(z, z)a(0)} 2 = ^^ + m + t 2 (^- + ±) + t ^n + --^ 1 (5.35) 

where fl = log(|) + 7, R oc r, and t, the scaling variable, is proportional to 
mr. 



We can use the zeroth order terms on both sides of |5.35| to fix R in terms 
of r. We find 

R = 2r (5.36) 
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At the next order, 6.351 reduces to 



in 



1 



R 1 / 2 4tt 



(5.37) 



were L(l) equals 



L(l) 



-A / d 2 «)(e i$ W/ 2 e i$ ( )e- $ W) 



(5.38) 



As with the first order term of the fermion correlator (ip^(z, z)ip + (0)), this 
integral has a logarithmic divergence. Regulating L as before by taking 
L(l) -> L(/3), we find 



L0) = ~\\zf /2 \z\ 2 ~ 2 ^ d 2 w\w-l 



\w\ 



(5.39) 



This integral has been done previously (see [4.27] ) with the result 



L0) = 2X\zf /2 \z\ 2 - 202 s0 2 /2)B(l - /3 2 /2, 1 - /3 2 /2)B(l - /3 2 /2,/3 2 - 1). 

(5.40) 

Taking p 1 to 1 via p 2 = 1 + e we find 



L (l + e ) = -27rm|^|^V a "V«| 2 ~ 2/5a ( ~ + 41og(2) 



(5.41) 



As before we do not expand \zf^ 2 /i^ 1 so as to preserve the correct scaling 
dimension. We then have 



L(l) =4vrm|2| 1 / 2 (log(^)+ 7 ). 



(5.42) 



Taking the scaling variable, £, to equal mr we see L{\) = -^rpiQ- Comparing 



this with |5.35| we see we have agreement. 

At second order, the perturbative expansion of the spin correlator gives 



RV 2 



1 



■M{1) 



(5.43) 
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where M(l) equals 
M(l) 



A 2 / ^^(e^e^W^e^e^) 



= X 2 \z\ s/I J d 2 wd 2 y\w-y\- 2 \w- lH^I- 1 ^ - lHyl. (5.44) 

We will show [5.43| is valid to the leading log term. 

M(l) is (quadratically) divergent. Again, as with the second order term 
of the fermion correlator, explicit inclusion of the bubble contribution leaves 
M(l) logarithmically divergent. But, as before, explicit inclusion of the 
bubbles is unnecessary. To regulate, we continue $. M{p) is then 

M0) = \ 2 \z\^ 2 ' 2 [ d 2 wd 2 y\w-y\-^ 2 \w-lf 2 \w\-^\y-l\^ 2 \yf 2 . (5.45) 



We have already evaluated this integral. Copying the result in [4.181 gives us 



M{(5) 



-4A 2 |z| 4 - 5/3 12 [J+ (s 2 {(3 2 /2)K_ + s([3 2 /2)s(3f3 2 /2)K + 
+J-(s 2 2 /2)K + -s 2 (f3 2 /2)K^ 



(5.46) 



where 
K ± 0) 

J±0) 



±B(1 ± P 2 /2, P 2 - l)r(l - P 2 W0 2 - 1)( y)0 2 - 1) x (5.47) 

3 F 2 (l±P 2 /2j 2 J 2 -lJ 2 ±P 2 /2,2,l) 
B(l ± P 2 /2, 2 - /3 2 )B(1 - P 2 , 1 ± /3 2 /2) x (5.48) 
3 F 2 (±/3 2 /2, 2 - /? 2 , 1 ± /3 2 /2, 3 - /3 2 ± /3 2 /2, 2 - /3 2 ± ^ 2 /2, 1). 



Setting p 2 = 1 + e we find 

1 



#±(e) 



T 2e 2 

1 r 



l + e(2- 7 -^(l±l/2)) + 0(e 2 ) 



1 + e(ip(l ± 1/2) + 7 =F 1 ± tt 2 /4) + 0(e 2 



(5.49) 
(5.50) 



where if) is the logarithmic derivative of the gamma-function, i.e. if>(x) 
T'(x)/T(x). The details of this calculation may be found in appendix E. 
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Putting everything together we find for M(l + e) 

M (l + e) = -A^-^^^-A-^ + O^e 2 )] (5.51) 



2tt 



2 



f 2 



z\- p2/2 \zm\ 2 (l -2elog(z/i) + 2e 2 log 2 



where in the last line we have substituted — m/i 1_/32 for A. The presence 
of the 1/e 2 term indicates, as expected, that the leading log term is log 2 . 
Substituting in t and _R 1//2 , subtracting the infinite terms, and dropping all 
but the log 2 terms, we find 

M(l) = -^2 (2vr 2 fi 2 + O(O)) . (5.52) 



Comparing with 5.35| we find the coefficients of the leading log's at second 



order agree. 

The ability to reproduce the non-trivial behavior of the Ising spin cor- 
relators provides us with a degree of confidence that our expression for the 
correlators away from (3 = 1 are correct. 

5.2 Gross-Neveu SU(2) 

We now go on to demonstrate that our methods of evaluating the sine-Gordon 
correlators reproduce the known behaviour of the Gross-Neveu model (sine- 
Gordon at p 2 = 2). Specifically we demonstrate that we are able to reproduce 
the known /^-functions for A and f3 governing the Kosterlitz-Thouless transi- 
tion to lowest order, as say calculated by Amit et al. |T3] and Boyanovsky 







As shown in section 3, the sine-Gordon theory develops UV singularities 
as p 2 — > 2. As with Ising at $ — 1, these singularities manifest themselves as 
poles in gamma functions. And similiarly to Ising, we regulate the correlators 
by continuing away from p 2 = 2 via p 2 = 2 — e. However unlike Ising, we 
handle these singularites, because they are UV, with a conventional renor- 
malisation. 



In order to facilitate comparison with Amit et al. [I2fl and Boyanovsky 



[O] we employ a sine-Gordon action equivalent to theirs: 

(5.1) 



S = -— I d 2 z 

47T 



d z <S>dz® + ^cos(/3$) 
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where the free propagator is 



($0)$(0)) = -log(4z2). (5.2) 

The difference in the propagators amounts to a difference in how its IR 
singularities are cured. 

To renormalise the theory we employ the following renormalisation pre- 
scription 

$o = Z*& R - 

# = Z^P R - (5.3) 
(cos(/3$)) = Z c (coa(fo)) R ; 

= Z c Z\Z$ \ R \ 

where the bare quantities are listed on the left hand side and the renormalized 
quantities on the right hand side. \r is understood to be dimensionless. This 
prescription, like Amit et al.'s and Boyanovsky's, involves a trivial renormal- 
isation of the (3 so as to ensure /5o$o = PrQr- However it differs from their's 
in that it introduces a wavefunction renormalisation of the cos(/3$). That 
this renormalisation is necessary is evident from considering the correlator 
(cos(/3$(^)) cos(/3$(0))). It possesses singularities at 0(A 2 ), the only way 
which to remove is through a wavefunction renormalisation. A consequence 
of this is that the /3-function for A has been computed incorrectly by both 



these authors. (We note that Amit et al.'s \12\ calculation of (3\ at higher 
orders was already known to be flawed because of his method of regulating 
the theory. See C. Lovelace |L4| for a detailed discussion.) However as we 



are interested in reproducing the /3-functions only at lowest order as a check 
on our methodology, we will not concern ourselves here with correcting this 
error. 

The wavefunction renormalisation for <3>, is found through calculating 
the correlator 

(d g $(z)d g $(0)). (5.4) 

Because this correlator is even in A, Zq, takes the form Zq> = Y^=o a 2nA 2n . 
Wavefunction renormalisation of cos(/?$) may be computed through the cor- 
relator (cos(/3$(z)) cos(/3$(0))). Knowing the wavefunction renormalisations 
of cos(/3<3?) and of $ allows the computation of the coupling constant renor- 
malisation, Z\, through the correlator (d z Q cos(/3$)). Because (d z Q cos(/3$)) 



36 



only has terms odd in A (and is finite at 0(A) ), Z\ takes the form 

oo 

Z x = ^(l + ]Ta 2n A 2 ™). (5.5) 

n=l 

H is an arbitrary mass parameter introduced to make the renormalised cou- 
pling dimensionless. Here we will not be interested in computing the higher 
order terms. So we have (rather trivially) 

Z x = ^(l + 0(A 2 )), (5.6) 

and we are left to compute Zq>. (We do not need to worry over Zq as it 
contributes to the renormalisation of A at 0(A 2 ), and so makes a contribution 
to A's /3-function only at 0(A 3 ).) 

To compute we employ a trick to calculate (d z $(z)d z &(0)). We can 
write 

(dMWMO)) = Jm^^^We^"))!^. (5.7) 

Thus we can use our results from Section 4. Doing so and taking the a — > 
limit we find 

(d z <S> (z)d z MQ)) = -{z)~ 2 - (5-8) 

Ar 



v 



16 • 2 2 WP% 



where L is 



L = -8np R s{p R ){p R - 1) — — ^— . (5.9) 

1 (3 — Pr) 

Writing /3 R — 2 — e, the singularities in the correlator become 

(0,$(z)0,$(O)> = -(^)- 2 + ^(^)-^(^) 2 -^ (1 + e(l/2 + 2 log(2))) . 

(5.10) 

Using minimal subtraction, we find Z$> to be 

^ = l-^ + 0(A 4 ). (5.11) 
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Knowing Z<$> and Z\ allows us to calculate the /3-functions for (3 and A: 



dX R dfi 
H^r~ = Mo 



djjL 



dpi 



-eX R + 0(X 3 



dZ, 



1/2 



dfj, 



64 



+ 0(A 4 ). 



(5.12) 
(5.13) 



Defining 5 = — e/2, we can recast the above in a form directly comparable 
with Amit et al.'s (12| and Boyanovsky's JO]] result 



25X R + 0{X S 



[3 R X* 



64 



+ 0(A 4 



^ + ^ + 0(A 4 ). 
32 32 v ; 



(5.14) 
(5.15) 



The leading terms of these /5-functions agree with their previous results. 
Thus we again have shown that this methodology is consistent with other 
methods. 



6 Discussion 

We have shown how the short distance expansion of some sine-Gordon corre- 
lation functions is well defined and reasonably tractable. The ultimate goal is 
to find some non-perturbative characterization of these correlation functions. 
We make two remarks. First of all, unlike the large distance expansion, the 
short distance expansions at and away from the free fermion point are of 
comparable complexity. Secondly, we have in no way utilized the integrabil- 
ity structure of the sine-Gordon theory, i.e. factorizable S-matrix, quantum 
inverse scattering method, etc. in developing the short distance expansion. 
It would be very interesting to understand the consequences of integrability 
in this context. 

The results presented herein leave two concrete avenues open for further 
exploration. The first arises from our knowledge of the differential equa- 
tion the correlators of the vertex operators satisfy at the free fermion point. 
Knowing the first set of terms in the perturbative expansion should allow, 
through the differential equation, the generation of higher order terms in the 
expansion. 
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The second is a correct evaluation of the /3-function at 0(A 3 ) for the cou- 
pling A in the SU(2) Gross-Neveu model. To make this evaluation two cal- 
culations would be necessary. The wavefunction renormalisation of cos(/3$) 
would need to be calculated through the correlator (cos($$(z)) cos(/3<&(0))). 
Knowing this then allows the calculation of Z\, a piece of the coupling con- 
stant renormalisation, at 0(A 2 ) through the correlator (d z <&(z) cos(/3$(0))). 
From there, f3\ is directly calculable. 
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A A Generalized Euler Integral 

In this appendix we evaluate the following generalized Euler integral, 

I = f 1 dw 1 w c t 1 (l-w 1 ) bl x (A.l) 
Jo 



dw 2 w a 2 2 {l - w 2 ) b2 {w 1 - w 2 ) ai2 x 

/■W2 

/ dw 3 w^ 3 (l — w 3 ) 3 (wi — w 3 ) ai3 (w 2 — w 3 ) 
Jo 



"23 ^ 



fWn-1 

/ dw n w a n n {l - w n ) bn (w 1 - w n ) ai " ■ ■ ■ (iu n _i - w n ) an -^. 
Jo 

Employing the following series of changes of variables, 

w 2 = w 2 wi, w 3 = w' 3 w 2 wi, ... , w n — w' n ■ • • W 2 Wi, 
the above integral reduces to 

/ = f'dw^^^^^il-w^ x (A.2) 
Jo 

f 1 dw.wf^^^^^il - w lW2 ) b2 (l - w 2 r 2 X 
Jo 

I' rf^3 EL3ai+E ^ < "" +n " 3 (l - W 1 W 2 W 3 ) b3 (l - W 2 W 3 ) a » X 

Jo 

(1 - w 3 ) a23 X 



\OC2n 

J n) 



f dw n w^(l- Wl ---w n ) bn (l-w 2 ---w n ) ain (l-w 3 ---w r , 
Jo 

(1 - W r ^ lWn ) an - 2 ' n {l - Wn) "' 1 '". 



Using the expansion, 



(1 - = £ L^V, (A.3) 



where (x)k = x{x + 1) • • • (x + k — 1) = F(x + k)/F(x), the above becomes 
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[ 1 dw 1 wF^ +E:<jaij+n ~ 1 (l-m) bl x (A.4) 
Jo 

y^tMhi / 1 ^"'^<--^ + ' , - 2 (i-, 2 )-(, 1 , 2 )^ x 

V ^W^ka / 1 ^ 3 ^3-^< !<3 ^ + -3 (1 _ w3r23 x 
1~ kl23 ] -k23 ] - JO 



fc 23 



^123^23! 

(wiw 2 w 3 ) hl23 (w 2 w 3 ) h23 X 



E (-&n)fci... n (-«ln)fc 2 ... n (-«2n)fc 3 ... n ' ' ' (-Qn-2,n)fc n _i,„ 
h \h \ h I 
fc 1 ... n "'l---n-" , 2---n- rtn,n— 1- 



x / dw n <" (1 - w^) ™-"- 1 {wi ■ ■ ■ w n ) kl - n ■ ■ ■ (w n w n - 1 ) kn ' n - 1 . 
Jo 

where the k m ... n arise from the expansion of (1 — w m ■ ■ -w n ). Defining the 
following notation 

n n 

<?m = J2 a *+ Yl <Xij + (n-m), (A. 5) 

i=m i^j>m—l 

k = sum of k's with index m (A. 6) 

m 

■■m---n + 

+ k 2 . 

■■m,m+l 

+ --- + k 2 . + 



km,m+l ^m,m+l,m+2 ' ' ' k m ... n7 

and 

(k)\ = k 12 \k 123 \---k 1 ... n \xk 23 \---k 2 ... n \x---xk n>n - 1 \, (A.7) 
I simplifies to 
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I = EIK- 6 i)*i-i U(-»ij)K + ,.,/(k)\x (A.8) 

{k} i=2 i<j-l 



f 1 dw 1 w a l 1+El \l-w 1 ) bl X 
Jo 

f 1 dw 2 w a 2 2+ ^ 2h (l-w 1 ) ai2 X 

JO 



f 1 dw n Wn n+T,nk (l ~ W 1 ) a - 1 ^ 

Jo 

Doing these integrals, we are able to put I into its final form: 

/ = B(6i + l,(7i + l)B(ai 2 + 1,(7 2 + 1) •• ■B(a n -i, n + l,<r n + 1) x 

n n 

EIK-kW..* II (-««Wi...,/(*) ! x 

{fc} i=2 «<j-l 

( <Tl + l) Eifc /( (Tl + 6 1 + 2) Eifc x (A.9) 
(^2 + l) S2fe /(^i + ai2 + 2) E2fe x 

(<r„ + l) En fc / ((7 n + a n _i,„ + 2)^ 2 fc . 

B Evaluation of K(i,j) 

In this appendix we evaluate the integrals K(i,j) introduced in Section 4.2. 
Using Appendix A we can write K(i,j) as 

K(i,j) = J* dw iW - 2+ P 2 (l - w t y' jT' dw jW ] 2+ P\l - Wjpiwi - Wj )-P 

= B( 7i + l,/3 2 -2)B(l-/3 2 ,/3 2 -l) x (B.l) 
y (-7,-) fc 2 ~ 2) fc 2 - l) fc 
fc =o A;!( 7l + /3 2 -l) fe (0) fc 

where B(x, y) = T(x)T(y)/T(x + y), (a)k = T(a + k)/T(a), and the sum J2k 
is the standard hypergeometric function, 

3 F 2 (- 7j , /3 2 - 2, /3 2 - 1, 7 , + P 2 - 1, 0, 1). (B.2) 
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However this form is somewhat problematic as a pole in 3 F 2 is multiplying a 

zero in B. To make sense of this expression we continue w~ 2+l3 to wj 2+l3 +e 
in the above integrand and take the limit e — > 0. With this continuation, 
K(i,j) equals 

K(i, j) = B( 7l + 1, ft - 2)B(1 -ft,ft-l + e)x (B.3) 
3 F 2 (-7„ /? 2 - 2, ft - 1, 7 , + /3 2 - 1, e, 1). (B.4) 

Using the identity 

ci m n f(0 = (R5) 



q(q + 1) • • • (a + w)/?(/3 + !)•••(/? + 71)7(7 + 1) ■ ■ ■ ( 7 + n) 

(n + l)!(J((J + l)...(<J + n) 
2 n+1 3 F 2 (a + n + 1, /3 + n + 1, 7 + n + 1, 5 + n + 1, 2, z), 



and then taking the limit e — > 0, we obtain 

K(i,j) = B(l + 7 ,,/3 2 - l)r(l - ftnP 2 - l)(-7i)C9 2 - l)x 

3 F 2 (l-7 J , / 9 2 -l,/3 2 ,/3 2 + 7j ,2,l), (B.6) 

where we have also used the relation xT(x) — T(x + 1). 



C Analytic Continuation of 3F2 

The Gaussian hypergeometric function 3 F 2 (q;, /3, 7, 5, 1) is, in general, a 
meromorphic function of the variables a, /3, 7, <5, and (. However its sum 
representation (as say given in Appendix B) is not necessarily convergent for 
all values of these variables. The convergence of the sum is determined by 
defining the "excess" parameter, s, by 

s = 5 + (-(a + (3 + j). (C.l) 

If s > the sum converges. However we will often encounter situations where 
s < and the sum representation is no good. To get around this problem 
analytic continutation formulae are needed. The only one that will prove to 
be necessary is an identity known as Dixon's theorem: 
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3 F 2 (a,f3, 1 ,5,(,l) = r( L { P T ^ C Fr( S \ S sF 2 (6 -a,(-a,s,s + (3,8 + ^,1). 
T(a)T(s + p)T(s + 7) 

(C.2) 

The convergence of the hypergeometric function of the right hand side is then 
determined by the positivity of a. Notice that because 3 F 2 is invariant under 
permutations of a, (3, and 7 and permutations of 5 and (, variations (144 in 
total - not all necessarily distinct) of this identity exist. 



D Taking $ 1 in (ip_(z, z)ip+(0)) 

We focus on the Taylor series expansions of the hypergeometric functions in 
K 2 (±(3). The expansion of the 3 F 2 function in K 2 (f3) is given by 



1 1,1 ' ' ' ' ' r(l + /3 2 )r(3-2/3 2 )(3-2/^ 

3 F 2 (2 - /3 2 , 1 - P 2 , 2 - /3 2 , 2, 4 - 2/3 2 , 1), (D.l) 



where we have used Dixon's theorem, C.2, and the identity, B.2. Evaluating 
further we have 

3F 2 (l + / 3 2 , / 3 2 -l,2-^ 2 ,2,l,l) = x 
V r(l + /3 2 )r(3-2/3 2 )(3-2/3 2 ) 

3 F 2 (2 - ^ 2 , 1 - ^ 2 , 2 - (3 2 , 2, 4 - 2/3 2 , 1) (D.2) 

= e + 0(e 2 ). 

The 3 F 2 function in K 2 {—(3) is simpler to evaluate: 

3 F 2 (l-/3 2 ,/3 2 -l,2-/3 2 ,2,l,l) = 3 F 2 (-e,e,l + e,2,l,l) 

= l + 0(e 2 ). (D.3) 



E Taking /3 1 in (a(z, z)a(0)) 2 

As in appendix D, we focus only on the hypergeometric functions that appear 
in the expressions that are being expanded about j3 — 1. Appearing in K±({3) 
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we have the function: 

Sk ± 0) = sF 2 (l ± /3 2 /2, /3 2 , /3 2 -lj 2 ± /3 2 /2, 2, 1). (E.l) 
To 0(e), S K± 0) is 

S K± 0) = 3 F 2 (1 ±1/2,1, e,l ±1/2,2,1) + 0(e 2 ) 

^ (2)**! 

r(2- e )r(i) 1 J 

= l±e±0(e 2 ), (E.2) 
where we have used the Gaussian summation formula 

^ (c) fc fc! r(c-a)r(c-6) 

In J ± we face the sum 

Sj ± 0)= 3 F 2 {±p/2,2-$ 2 ,l±p/2,3-$ 2 ±p/2,2-p±p/2,l). (E.4) 

Using the analytic continuation formula in appendix C, Sj ± (f3) can be trans- 
formed into 

r r(3-^±^/2)r(2-^ 2 ) . 

5j±(/5) " r(i±^/2)r(4- 2 ^) SjM (E - 5) 

where Sj ± {(3) is 

S J± = 3^(1 - P 2 , 2 - /? 2 , 2 - ^ 2 , 2 - (3 2 ± ^ 2 /2, 4 - 2^ 2 , 1). (E.6) 
To 0(e), S J± 0) is 

Sj ± (e) = 3 i 71 2 (-e,l,l,l±l/2,2,l) + 0(e 2 ) (E.7) 

(l) fc _i(l)fc(l)* 



jfefc!(l±l/2) fc (2) fc - 
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The evaluations of the two sums in |E.7| are simlar. So we will only be explicit 
in the evaluation of Sj + (e). 
We may write 



Sjje) 



(i)*-i(i)fc 



i (3/2) fc (2), 



(E.8) 



1 — e lim 



(1 + S) k - 2 (1 + S) k . 



k=2 



(3/2) fc _i(2) fc . 



i 



e T(5-l)T(5) ^ (5-l) k (5) k 

2^0 r*(i + S) t 2 (i/2) k k\ ' 



Now using the Gaussian summation formula we find 



SjAe) = l-~lim 



T(S-1)T(S) 



2^o T 2 (l + 5) 
Evaluating this limit we obtain 



T(3/2 - 25)r(l/2) 
r(3/2 - 5)r(l/2 - 5) 



1-26(5-1) 



,7T 



5 J+ (6) = l + e ( T -3)+0(6 2 ). 



(E.9) 
(E.10) 



Similarly we find for S j_ (e) 



,7T 



^_(e) = l-e( y + 2) + 0(e 2 ). 



(E.ll) 
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